Topological phase transitions in ultra-cold Fermi superfluids: 
the evolution from BCS to BEC under artificial spin-orbit fields 
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We discuss topological phase transitions in ultra-cold Fermi superfluids induced by interactions 
and artificial spin orbit fields. We construct the phase diagram for population imbalanced systems 
at zero and finite temperatures, and analyze spectroscopic and thermodynamic properties to charac- 
terize various phase transitions. For balanced systems, the evolution from BCS to BEC superfiuids 
in the presence of spin-orbit effects is only a crossover as the system remains fully gapped, even 
though a triplet component of the order parameter emerges. However, for imbalanced populations, 
spin-orbit fields induce a triplet component in the order parameter that produces nodes in the 
quasiparticle excitation spectrum leading to bulk topological phase transitions of the Lifshitz type. 
Additionally a fully gapped phase exists, where a crossover from indirect to direct gap occurs, but a 
topological transition to a gapped phase possessing Majorana fermions edge states does not occur. 

PACS numbers: 03.75.Ss, 67.85.Lm, 67.85.-d 
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Ultra-cold Fermi atoms are one of the most interesting 
physical systems of the last decade, as they have served as 
quantum simulators of crossover phenomena and phase 
transitions encountered in several areas of physics. Due 
to their tunable interactions, atoms like ^Li and have 
been used to study the crossover from BCS to BEC super- 
fluidity, to simulate superfluidity in neutron stars, and to 
investigate unitary interactions which are of great inter- 
est in nuclear physics. Furthermore, the ability to con- 
trol the internal spin state of the atoms by using radio- 
frequencies (RF) enabled the studies of quantum and 
classical phase transitions as a function of interactions 
and population imbalance. These tools have permitted 
the study of crossover phenomena and phase transitions, 
and have validated the symmetry based classification of 
phase transitions put forth by Landau over the thermo- 
dynamic classification proposed earlier by Ehrenfest. 

Very recently a new tool for the toolbox was created: 
artificial spin-orbit coupling has been produced in neu- 
tral bosonic systems [l| where the strength of the cou- 
pling is controlled optically from weak to strong. The 
same technique can be applied to ultracold fermions [3| 
and should allow for the exploration of superfluidity as 
a function of interactions and fictitious spin-orbit cou- 
pling 0. This possibility created enormous theoretical 
interest recently 0^ , which was focused on zero temper- 
ature (ground state) properties. Considering possible ex- 
periments with fermionic atoms such as ^Li, '*°K, we dis- 
cuss here topological phase transitions at zero and finite 
temperatures for imbalanced fermions during the evolu- 
tion from BCS to BEC supcrfiuidity in three dimensions 
and in the presence of controllable spin-orbit couplings. 
Even though the symmetry of the order parameter does 
not change through topological phase transitions, violat- 
ing Landau's symmetry-based classification, clear signa- 
tures emerge in spectroscopic and thermodynamic prop- 
erties validating Ehrenfest's ideas, which combined with 



changes in topological (Hopf) invariants produce a finer 
classification scheme of phase transitions. 

Hamiltonian: We start with the Hamiltonian density 

H{r) = Hoir)+Hi{r), (1) 

where the single-particle term is simply 

■Ho(r) = (r) kaSap - hi{v)a^^ap V'/3(r)- (2) 

Here, Ka ~ — V^/(2mQ) — /.ic is the kinetic energy in ref- 
erence to the chemical potential , and hi (r) is the spin- 
orbit field along the i-direction [a =t, J,, i = a;, y, z). The 
interaction term is 'Hi{v) = —g^].{r)'il?^{r)tl>^{r)tl>^{r), 
where g is a, contact interaction, and we set h = = 1- 
Effective Action: The partition function at tempera- 
ture T is Z = / VIiIj, exp {-S[tp, with action 



tl - 



drdr 



(3) 

Using the standard Hubbard-Stratanovich transfor- 
mation that introduces the pairing field A(r,T) = 
g{'ipl{r^ r)'0^(r, r)) and integrating over the fcrmion vari- 
ables lead to the effective action 

|A(r,r)p 



Sctf 



drdr 



T , , M 

in det — 

2V T 



K+S{r - r') 



where K+ = (Kf + K^)/2. The matrix M is 



M 



( dr+ A> 

At - K-^ 



-A \ 







V 



-At 







(4) 



h±_ dr - R'l J 

where h± = /ij. — ihy corresponds to the transverse com- 
ponent of the spin-orbit field, to the parallel com- 
ponent with respect to the quantization axis z, = 
- hz, and ATj, = A'j, + h^. 
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Saddle Point Approximation: To proceed, we use the 
saddle point approximation A(r, r) = Ag + r](r,T), and 
write M = Mq + IsAp. The matrix Mg is obtained via 
the saddle point A(r,T) Aq which takes M — ^ Mq, 
and the fluctuation matrix Mp = M — Mq depends only 
on 77(r,T) and its Hermitian conjugate. Thus, we write 
the effective action as 5'off = Sq + Sp. The first term is 



Sn 



V\^o\ 
T a 



T 



k 



in momentum- frequency coordinates (k, where 
bJn = (2?! + 1)ttT. Here, Ej{k) are the eigenvalues of 



H 



— 



-/^l(k) A';(k) 



V 



-A 



which describes 
tary excitations 



a; 



the 





-A>(-k) 
/^±(-k) 



vfT = {^t(k),V;j(k),Vt(-k),Vi(-k)} 
ation action is 

>(r,r)p 



-Ao 


/il(-k) 

(5) 

Hamiltonian of elemen- 
tlie four-dimensional basis 
The fluctu- 
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In dot (l-hMo^Mp) 



The spin-orbit field is h^(k) = hfl(k) -I- hpfk ), w here 
hi<'(k) = vpiX—ky'k + kxy) is of Rashba-type Il0| and 
h.r)l\i) = vd {kySi. + k^y) is of Dresselhaus-type has 

magnitude /i^(k) — ^ {vd ~ vnf fc^ _|_ ^ 

For Rashba-only (RO) (ud = 0) and for equal Rashba- 

Dresselhaus (ERD) couplings {vr = vd = v/2), the 

transverse fields arc /i±(k) — VR^k'^ + k^ {vr > 0) and 

h±{'k) ~ v\kx\ {v > 0), respectively. 

Order parameter and number equations: The thermo- 
dynamic potential is 51 = fig + f^F, where 

no = - I ^ In {1 + exp [-E,ik)/T] } + Y.K+, 

with K+ = K^{-k) + Kii-k) /2 is the saddle 

point contribution and ^Ip = —TlnZp, with Zp = 
J D[r7, 77] exp [— S'f(^, 77)] is the fluctuation contribution. 
The order parameter is determined via the minimization 
of ilg with respect to |Agp leading to 



^ = -lY.^p[E,ik)] 



dE,{k) 



(6) 



where np [Ej{k)] = l/(exp [£'j(k)/T] + 1) is the Fermi 
function for energy Ej{k). We replace the contact in- 
teraction g by the scattering length Os through the rela- 
tion 1/g = -m+/(47ra,) + il/V)J2k [l/(2ek,+)] , where 



m+ = 2mim^/{mi + nif) is twice the reduced mass, 
Ek.Q = fc^/(2mc() are the kinetic energies, and ek.+ = 
[^k.t + £k.i] /2. The number of particles for each spin 
state is ~ —dfl/dna, which is written as 



No, = No,,0 + A^a,F, 

where the saddle point contribution is 



(7) 



dflo 



, „ dE,(k) 



and the fluctuation contribution is Na_p = —dflp/d^a- 
The self-consistent relations shown in Eqs. © and ([7]) 
are unphysical for atoms of unequal masses, as it is 
not possible to have a spin-orbit fleld that converts 
atom A into atom B, as this violates barionic num- 
ber conservation. Therefore, we particularize our dis- 
cussion to balanced and imbalanced systems of equal 
masses, where the explicit eigenvalues of the matrix 

Hg are Ei{k) = ^Jy+ +Y^+ 2^Y+Y^ - |Ao|2|^, 
which is always positive deflnite, and E2{k) = 



Y+ + Y-- 2y/Y+Y- - |Ag|2|/ij_|2, which can have ze- 
ros in momentum space, £'3(k) = — £'2(k) and £'4(k) ~ 
-Ei{k). Here, Y+ = Kl(k) + |Ag|2 > 0, r_ = Al(k) + 
/ij^(k)p > 0, with K± = {K^±K^)/2. For equal masses, 
= Ek ~ with = k^ /2m being the single particle 
dispersion and /i+ = {^f + ^i)/2 being the average chem- 
ical potential. Furthermore, = (p,^ — ni)/2 ~ fi^ 
plays the role of a Zeeman field that causes population 



imbalance P 



N^/N+, where 7V_ 



- Ni and 



A^_i_ = N-f + N^. We define momentum, energy and veloc- 
ity scales through the total particle density = n-^ + ui 
or fc|,_|./(37r2) = fc|,^/(67r2) + kl^/{6Tr^). This choice 
leads to the Fermi momentum kp_^^ = {kp^ + kp^)/2, 
Fermi energy ep^ ~ fc|,^/2TO and to the Fermi velocity 
vf+ = kpj^jm. 

Phase Diagram: The solutions of Eqs. ^ and ([7]) sub- 
ject to thermodynamic stability conditions (positivcness 
of compressibility matrix Rap = {dNa/dfip)T.v-, isother- 
mal compressibility kt and volumetric specific heat Cy) 
lead to the T = phase diagram of population imbal- 
ance P versus interaction parameter l/(kp^as), which 
is shown in Fig. [T] for various values of ERD spin-orbit 
coupling v/vpj^. The phase diagram for v/vp^ = 
shows a large region where a uniform superfiuid US phase 
is unstable, we call this region NU for possible non- 
uniform phases including phase separation between su- 
perfiuid and non-superfiuid components or modulated 
superfiuid phases. A general tendency of the phase di- 
agram with increasing spin-orbit coupling v/vp^ is the 
stabilization of superfiuid phases and the shrinkage of 
unstable regions. This stabilization is largely due to the 
emergence of a triplet component in the order parame- 
ter for superfluidity, which circumvents the pair breaking 
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FIG. 1: Phase diagrams of P versus l/(kF+as) for ERD spin- 
orbit couplings: a) v/vf+ = 0; b) v/vf+ ~ 0.28; c) v/vf+ ~ 
0.57; d) v/vf+ = 0.85. Solid lines are phase boundaries, 
the dashed line corresponds to ^-i- = 0, and the dotted line 
corresponds to a crossover. 



abilities of Zeeman fields /i_ . Many multi-critical points 
and phases exist, including a normal (N) and three uni- 
form superfluid regions: the phase US-0 corresponds to 
a fully gapped superfluid with no (zero) line of nodes, 
the phases US-1 and US-2 correspond to gapless super- 
fluids with one-line or two-lines of nodes, respectively. 
Within the gapped US-0 bulk phase, in our three dimen- 
sional system, a topological transition due to emergence 
of chiral edge states (Majorana fermions) does not occur, 
instead we have a crossover from indirectly to directly 
gapped US-0. 

Spectroscopic Properties: For the excitation energies 
i?i(k) and i?2(k), only -E2(k) can have zeros. In the 
ERD case, where /i±(k) = v\kx\, zeros of -E2(k) can oc- 
cur for y+ = Y- with fc^; = 0. This leads to the following 
cases: (a) two possible lines (rings) of nodes at (fc^ + 

kl)/{2m) 



Ai4 



|Aop for the outer ring, and 



{ky + k^)/{2m) = /i+ - y yu^ - |Aop for the inner ring, 

when /i^ — |Aop > 0; (b) doubly-degenerate line of nodes 
{ky+kl)/{2m) ~ /i+ for ^+ > 0, doubly-degenerate point 
nodes for /i+ = 0, or no-line of nodes for < 0, when 
— |AoP ~ 0; (c) no line of nodes when fit — |Aop < 0. 
The study of case (a) can be refined into cases (a2). (al) 
and (aO). In case (a2), two rings indeed exist provided 

that > \/ IJ-t — |Ao|^. However, the inner ring dis- 



appears when = y A'- ^ l^ol^- In case (al), there 
is only one ring when < ^ — |Ao|^. In case 

(aO), the outer ring disappears at ^+ = — J fit — |AoP, 



and for fi+ < — y/^^ — |Ao|^ no rings exist. Thus, 
the US-2/US-1 boundary is determined by the condition 
fi+ — fit ~ |Ao|^; the US-O/US-2 boundary is deter- 
mined by the Clogston-like condition \fi-\ = |Ao| when 
fi+ > 0, where the gapped US-O phase (singlet-rich) dis- 
appears leading to the gapless US-2 phase (triplet- rich) ; 
and the US-O/US-1 phase boundary is determined by 

fi+ = — ^ fit — I Aq I ^ . Furthermore, a crossover line be- 
tween an indirectly gapped and a directly gapped US- 
phase occurs for fi+ < and = |Ao|. Lastly, 

some important multi-critical points arise at the inter- 
sections of phase boundaries. First the point ^+ = and 
= |Ao| corresponds to a tri-critical point for phases 
US-0, US-1, and US-2. Second, the point |Ao| = and 
/i+ = corresponds to a tri-critical point for phases 
N. US-1 and US-2. Representative excitation spectra are 
shown in Fig. [21 Near the zeros of -E2(k), quasi-particles 
have linear dispersion and behave as Dirac fermions. The 
change in nodal structures is associated with bulk topo- 
logical transitions of the Lifshitz class as noted for p- 
wave [l^ and d-wave [l^l superfluids The loss of nodal re- 
gions correspond to annihilation of Dirac quasi-particles 
with opposite momenta, which lead to significant changes 
in spectroscopic properties such as momentum distribu- 



tions [13[ nk , 



l-E,MEjmdE,ik)/dfi^ /2 



o' 



o 




FIG. 2: Excitation energy E2(k)/eF+ versus k±/kF+, where 
k± = ^/kf+k?, for v/vF+ = 0.57: a) US-2 phase with P = 

0. 5 and l/{kF+as) = -0.92; b) US-1 phase with P = 0.5 and 
l/{kF+as) = 2.0; c) US-0 phase (indirect gap) with P = 0.06 
and l/(fcF-t-as) = —0.5; and d) US-0 phase (direct gap) with 
P = 0.06 and l/(fcF-f a,) = 0.5. 

Topological Order: The superfluid phases US-0, US- 

1, US-2 are characterized by different effective actions 
S'cff which depend explicitly on the matrix M(icj,k) = 

-Ho(k)]" • Setting iu) ^ Lo and using algebraic 
topology [IJI, we construct the topological invariant 



m 
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which in the gapped US-0 phase is to = 0, in the gapless 
US-1 phase is m = 1, and in the gapless US-2 phase is 
TO = 2, such that to counts the number of rings in each 
phase. The integral above has a hyper-surface measure 
dS^ and a domain 2? that encloses the region of zeros of 
Lu = EjCk) = 0. Here /x, A,7 run from to 3, and 
has components feg = ki = k^, fc2 = ky, and k^ = k^. 

Thermodynamic Properties: To characterize topolog- 
ical phase transitions thermodynamically, we show in 
Fig. |3] the compressibility matrix function of 

l/ikp+as) for P = 0.06,0.14,0.5,0.8 at ERD spin-orbit 
coupling v/vf+ = 0.57. Notice that singular behavior 
occurs at phase boundaries. For instance, in Fig. [3^, K^p 
has cusps at the phase boundaries US-2/US-0 and US- 
0/US-l and are strictly positive. However, in Fig. [3t, 
RajS has discontinuities at the N/US-2 boundary, diver- 
gences at the US-2/NU and NU/US-1 boundaries and 
the diagonal elements Raa are negative in the NU region. 



superfluids and described multi-critical points. In addi- 
tion, we have found that spin-orbit effects tend to stabi- 
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FIG. 3: The T = compressibility matrix ep+hiai) / N+ as a 
function of 1/ [kp+as) for v/vf+ = 0.57 at population imbal- 
ances a) P = 0.06, b) P = 0.14, c) P = 0.50, d) P = 0.80. 

Finite Temperatures: As experiments are not per- 
formed at T = 0, we show in Fig. U] the P versus T phase 
diagram at unitarity 1/ {kF+ag) = for various ERD 
spin-orbit coupling: a) v/vf+ = 0; b) v/vf+ = 0.28; c) 
v/vf+ = 0.57; and d) v/vf+ = 0.85. Phase boundaries 
are determined by singular behavior of In Fig.[4l the 
reduction of the NU region with increasing v / vf+ is due 
to the increasing importance of the triplet component of 
the order parameter, which stabilizes phases US-1 and 
US-2 at intermediate and large P. 

Conclusions: We have investigated spectroscopic and 
thermodynamic properties in the BCS-to-BEC evolution 
at zero and finite temperatures, including spin-orbit ef- 
fects. We have identified bulk topological phase transi- 
tions of the Lifshitz class between gapped and gapless 



FIG. 4: Phase diagrams of P versus reduced temperature 
T/Tc{v) at unitarity l/{kF+as) = for: a) v/vf+ = 0; b) 
v/vF+ = 0.28; c) v/vF+ = 0.57; d) v/vf+ = 0.85. 



lize uniform superfluid phases against phase separation 
or non-uniform superfluid phases as population imbal- 
ance or interactions are changed, due to the creation of 
a triplet component in the order parameter that circum- 
vents pair breaking effects due to Zceman fields. 
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